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Abstract 

We present certain mathematical aspects of an information method 
which was formulated in an attempt to investigate diffusion phenom- 
ena. We imagine a regular dynamical hamiltonian systems under the 
random perturbation of thermal (molecular) noise and chaotic motion. 
The random effect is taken into account via the uncertainty of irregular 
dynamic process produced in this way. This uncertainty due to differ- 
ent paths between two phase points is measured by a path information 
which is maximized in connection with the action defined originally for 
the unperturbed regular hamiltonian systems. The obtained transi- 
tion probability depends exponentially on this action. The usefulness 
of this information method has been demonstrated by the derivation 
of diffusion laws without the usual assumptions. In this work, some 
essential mathematical aspects of this irregular dynamics is reviewed. 
It is emphasized that the classical action principle for single least ac- 
tion path is no more valid and the formalism of classical mechanics for 
regular hamiltonian systems is no more exact for irregular hamiltonian 
dynamics. There is violation of the fundamental laws of mechanics by 
randomly perturbed hamiltonian systems. However, the action prin- 
ciple is always present for the ensemble of paths through the average 
action. This average action principle leads to a formalism of stochas- 
tic mechanics in which, in spite of the violation of fundamental laws, 
the mathematical form of classical mechanics can be recovered by a 
consideration of the statistical averaging of the dynamics. 

PACS numbers : 02.50.-r (Stochastic processes); 66.10.Cb (Diffusion); 
45.20.-d (classical mechanics) 
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1 Introduction 



The dynamics of a Hamiltonian system can be roughly classified into two 
categories: regular and irregular. The mathematics of the regular dynamics 
can be perfectly formulated on the basis of the least action principle in clas- 
sical mechanics. But the stochastic character of irregular dynamics makes it 
much more complicated to be formulated. Description of random dynamic 
process with Newtonian type equations is possible if one introduces random 
forces like in Langevin equation [1]. Other statistical method by informational 
consideration without introducing the details of the dynamic randomness are 
proved to be very useful. The description of the stochastic behavior of chaotic 
systems using Kolmogorov-Sinai entropy is an example[2]. 

The method we address in this paper is originally intended to study dif- 
fusion process by a informational consideration. Diffusion is an irregular 
dynamic process in which components of a mixture are transported around 
the mixture by means of random molecular motion. Over 200 years ago, 
Berthalot postulated[3] that the flow of mass by diffusion across a plane, 
was proportional to the concentration gradient of the diffusant across that 
plane. About 50 years later, Fick introducedfl] two differential equations 
that quantified the above statement for the case of transport through thin 
membranes. Fick's First Law states that the flux J of a component of con- 
centration n across a membrane is proportional to the concentration gradient 
in the membrane: 

= (D 

where x is the position variable for one dimensional systems. Fick's Second 
Law states that the rate of time change of concentration of diffusant at a 
point is proportional to the rate of spacial change of concentration gradient 
at that point within the mixture 



On d 
dt dx 



dn(x) 
dx 



(2) 



If D is constant everywhere in the mixture, the above equation becomes 

dn _ p ) d 2 n(x) 
dt ~ U dx 2 ■ 

Above normal diffusion laws are very precisely tested in experiments in 
most of solids, liquids and gases and are widely studied in nonequilibrium 
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thermostatistics together with the Fokker-Planck equation of diffusion prob- 
ability, the Fourier law of heat conduction and the Ohm's law of electrical 
charge conduction. Many efforts to derive theoretically the above diffu- 
sion laws were concentrated on special models of solids in which particles 
are transported. Other phenomenological derivations are also possible if 
one supposes Brownian motion and Markovian process [1], or Kolmogorov 
conditions [5]. Recently, an attempt is made to derive these laws with a 
different method based on information theory. The method consists in con- 
sidering hamiltonian systems under the perturbation of thermal noise and 
chaotic instability. The dynamic uncertainty of this perturbed dynamic pro- 
cess is measured by a path information connected with different paths. In 
order to derive the transition probability distribution of the irregular pro- 
cess, the path information is maximized in connection with, faute de mieux, 
the action defined with the unperturbed Hamiltonian. The diffusion laws is 
simply the differential equation of the transition probability distribution in 
exponential of action. Although we can hardly talk about the exact action 
of a system subject to random forces and that the action here is only the 
unperturbed action of the system, the least action principle is still present 
in this approach through the average action. This description of irregular 
dynamics underlies a probabilistic formalism of mechanics which is logically 
different from the conventional one based on the action principle for single 
path of the unperturbed hamiltonian systems. In what follows, we discuss 
some mathematical aspects of this stochastic mechanics. 



2 Regular dynamics of hamiltonian systems 

A hamiltonian system is a mechanical system which satisfies the following 
Hamiltonian equations [2] : 



x\ = — — and p\ = — — — , with i = 1,2, ...n (3) 
dpi dxi 

where n is the number of degrees of freedom, X{ is the coordinates, Pi = mx\ 
the momenta and H is the hamiltonian of the system. The solutions of these 
equations are geodesies in phase space Y of 2n dimension, i.e., the paths 
between a initial point a and a final point b whose action = J a L^{t)dt 
has a stationary (least action principle). From this action principle, we can 
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derive Euler-Lagrange equations [6] 



d dL k (t) _ dL k (t) = Q 
dt dx dx 

where L k = T — V (we suppose n = 1 from now on) is the Lagrangian 
of the system along a path k. If we consider the Legendre transformation 
H = Px — Lk, Eqs.(3) can be easily derived from Eq.(4)[6]. This formalism of 
mechanics underlies a completely deterministic character of dynamic process: 
if the initial conditions are given, there is no uncertainty, no information, no 
probability distribution associated with the future evolution of the systems. 
The geodesies are determined for any time by the above equations plus initial 
conditions. 



3 A path information 

Now we imagine a hamiltonian system under the perturbation of thermal 
noise and chaotic instability. This perturbation can be either internal or 
external to the system. A hamiltonian system containing a large number 
of particles in random motion is an example of internal perturbation. The 
dynamic instability of chaotic systems is another example of internal per- 
turbation leading to stochastic deviation from regular behavior. The motion 
of a Brownian particle is an example of external perturbation due to the 
random molecular motion around the particle. Under these perturbations, 
the geodesies of regular dynamics will be deformed and fluctuate in such a 
stochastic way that following exactly the evolution of each mechanical quan- 
tity (e.g., action) is inconceivable. We call this dynamic process irregular 
since there is uncertainty about the future evolution even if the initial con- 
dition is given. We will consider two dynamic uncertainties: 

1. Between any two phase space points a and b, there may be different pos- 
sible paths (labelled by /c=l,2,...w) each having a probability pk(b\a) to 
be followed by the system. This is the uncertainty considered by Feyn- 
man in his formulation of quantum mechanics[7]. Here we introduce it 
within randomly perturbed classical mechanics. 

2. There are different possible paths leaving the point a and leading to 
different final points b, each having a probability pk(x \ a) to be followed 
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by the system, where x is the position of arbitrary b. This uncer- 
tainty is the basic consideration for the definition of Kolmogorov- Sinai 
entropy[2] in the description of chaotic systems. 

In what follows, we first discuss the transition probability pk(b\a) between 
two fixed phase points (or cells of a given partition of the phase space): the 
cell a in the initial phase volume A and a cell b in the final phase volume 
B via a path k (k = 1, 2, ...w). We suppose that the uncertainty concerning 
the choice of paths by the systems between two points is measured with the 
following path information 

w 

hb = Pk{b\a) In p k (b\a). (5) 
fc=i 

We have the following normalization 

w 

5>*(6|a) = l. (6) 

k=l 

4 Maximum path information 

An essential assumption of our informational method is that each path is 
characterized by its action A ab {k) defined as if there was no random pertur- 
bation. This action is given by 

A ab (k) = f L k (t)dt (7) 

Jt ab (k) 

where Lfc(t) = E — U is the Lagrangian of the system at time t along the 
path k. The average action is given by 

A ab = Y / Pk(b\a)A ab (k). (8) 

k=l 

For irregular dynamic process, the least action principle does not apply since, 
firstly, there is no exact calculation of real action due to dynamic randomness, 
and secondly, the action we defined here has no stationary 5A ab {k) = 
for all the possible paths. We use the information method which consists 
in saying that the stable probability distribution of paths must correspond 
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to a stationary uncertainty or path information I ab under the constraint 
associated with the average action A ab . This means the following operation: 



w w 



S[I ab + aJ2Pk(b\a) - r]^p k {b\a)A ab {k)] = 



(9) 



k=l k=l 



leading to 



Pk(b\a) = —exp[-r)A ab (k)], 



(10) 



where the partition function 



(11) 



k 



The physical meaning of the multiplier rj will be discussed below. 

It is proved that [8] the distribution Eq.(10) is stable with respect to the 
fluctuation of action. It is also proved that Eq.(10) is a least (stationary) 
action distribution, i.e., the most probable paths are just the paths of least 
action. Indeed, we have 5pk(b\a) = —r]p k (b\a)5A ab (k) = which means 
5A ab (k) = and yields Euler-Lagrange equation and Hamiltonian equations 
given in Eqs.(3) and (4). Consequently, the most probable paths are just the 
geodesies of the unperturbed hamiltonian dynamics. The other paths do not 
satisfy Eqs.(3) and (4). In general, the paths have neither SA ab (k) = nor 
SA ab = 0. Nevertheless, the average action A ab defined over all these paths 
does have a stationary associated with the stationary of the uncertainty I ab 



Here we considered J2\"=i $Pk(b\a) = 0. Eq.(12) implies that, although the 
least action principle of classical mechanics cannot apply when the dynamics 
is perturbed by random and instable noise, the maximum path information 
introduced above underlies in fact the same physics in which the action prin- 
ciple is present as a average effect in association with the stationary dynamic 
uncertainty. Eq.(12) implies a probabilistic version of classical mechanics 
and will be used below to derive the averaged Euler-lagrange equation for 
irregular dynamics. 



in Eq.(9), i.e., 



- r)SA ab + SI ab 



0. 



(12) 
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5 A calculation of least action distribution 



In what follows, the unperturbed action is analyzed for different paths with 
the Euler method. Let us look at a system of mass m moving along a given 
path k from a point a to a point b. The path is cut into N infinitesimally 
small segments each having a spatial length Axi = Xi — Xi-i with i = 1...N 
(xq = x a and xn = Xb). t — U — U_\ is the time interval spent by the system 
on every segment. The Lagrangian on the segment % is given by 



X{, t) 



m(xi — Xj_i) 2 / dU" 
2(i i -t,_ 1 ) 2 \dx~ / 



(Xi 



U(x 



i-lj 



(13) 



where the first term on the right hand side is the kinetic energy of the parti- 
cle, the second and the third terms are the average potential energy on the 
segment i. The action of segment % is given by 



where Fa 



A, 



dx J , 



m(Axi 
2T 



^ Axi 



U(x 



i-l 



(14) 



Pk,i 



is the force on the segment i. According to Eq.(lO), the 

(15) 



transition probability pk t i from Xi-\ to Xi on the path k is given by 

1 



a exp 



-V 



—Ax- + Fi-Ax t 
It 1 2 



where Zj is calculated as follows 



dxi exp 
i 8m 



-V 



—Ax 2 + Fi-Ax, 
2t 1 2 



(16) 



exp 



'27rf 

mr\ 



The potential energy of the point x,_i disappears in the expression of p^, 
because it does not depend on X{. 
The total action is given by 



N 



1=1 



i=l L 



m(Axi] 
2t 



Fi-Axi-Uix^t 



(17) 
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According to Eq.(lO), the transition probability from a to b via the path k 
is the following: 



N 



Pk(b\a) 



exp -7] 



i=i 



2~T 



Fi-Axi 
2 



(18) 



JV 



p(b\a) 1 Y[p k , 



i=i 



where 

w / 

z = Y, ex P\-vY, 



k=l 

fOO 

oo 



N ^(AXi) 2 t A 

+ Fi- Ax, 



i=l 



2t 



(19) 



J k/ 



j j j ( ^ \m{xi - 

dx\dx2--.dxN-i exp —77 7 y 



2f 



+ -Fi-(xi - Xi-i) 



exp 



,2^ 
i 8m 



F2nt\ 



N 



I p(6|a) = Zft>(6|a) 



and 



p{b\a) = exp 
x exp 



>Mh - t a y 



5m 



777,77 



2ir(t b - t a 



(20) 



V 



m ( x b ~ Kg? F h^ta A 

2(4 "ta) 1 2 



Remember that in the above calculation, the point Xq = x a and the final 
point xn = Xb are fixed. 

Now in order to see the behavior of transition probability with respect 
to final point, we have to relax x\, = x and let it vary arbitrarily as other 
intermediate points. This implies we take into account the second uncertainty 
due to chaos mentioned in the introduction. The corresponding transition 
probability pk{x\a) from a to arbitrary x via the path k has been derived 
with the maximum path information combined with action[9]. Here we only 
introduce it in the following way: 



Pk{x\a) 



p(b\a)pk(b\a) 

N 

Y[pk,i/i-i, 

i=l 



(21) 
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normalized by 

w „ 

J2J2Pk( x b\ a ) = / dx 1 dx 2 ...dx N -idxp k (x\a) = 1. (22) 

b k=l J 

The uncertainty or path information associated with this transition proba- 
bility between a and arbitrary b is given by 

w 

I a = -^2^2pk(x\a)\np k (x\a) (23) 

b k=l 

w 

= -J2p( b \ a ) ln P( b \ a ) ~ J2p( b \ a ) J2Pk( b \ a ) ln Pk( b \a) 
b b k=l 

= h a b + (lab) 

where h a b = — J2bP( b \ a ) hxp(b\a) is the uncertainty of the transition from a 
to an arbitrary b via whichever path, and (I ao ) is the average of I a b over all 
the possible point b in the final phase volume B. 

The average action defined with pk(x\a) over all the possible paths and 
final points is given by 

w 

(A ab ) = -J2T,Pk( x \ a ) A *b( k ) ( 24 ) 

b k=l 

w 

= J2p( b \ a )T,Pk( b \a)Aab(k) 

b k=l 

= Y.p( b \ a ) A *b- 

b 

Remember that p(b\a) does not depend on the paths. So if A a t, has a sta- 
tionary, (A a b) must have a stationary as well. 

6 Fokker-Planck equation 

A derivation of the diffusion laws from Eq.(15) and Eq.(21) is given in [10]. 
In what follows, we only illustrate this straightforward method with Fokker- 
Planck equation which describes the time evolution of transition probability. 
This equation can be derived if we suppose that the diffusion particles follow 
Brownian motion and Markovian process [1], or Kolmogorov conditions [5]. 
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From Eq.(15), a simple calculation of the derivatives d ^ p ^ and Qffi 
yields 

ggj = r d(F^) | 1 (9 2 ^ 

This is the Fokker-Planck equation, where r is the mean free time supposed 
to be the time interval t of the particle on each segment of its path. In view 
of the Eq.(21), it is easy to show that this equation is also satisfied by pk(x\a) 
if Xi is replaced by x, the final position. Other diffusion laws can be easily 
obtained from the above equation. We think that this derivation is a proof of 
the usefulness of the present informational approach to irregular dynamics. 
As discussed in [10], Eq.(25) yields in general: 

dn t d{nF) ^ 1 d 2 n 
dt m dx 2mr] dx 2 

where n is the particle density at point x and time t. In order to clarify the 
physical meaning of 77, we consider here a gas in equilibrium (hence |j| = 0) 
in a constant force F with the distribution n(x) = n(xo)exp[F(x — xo) / ksT] 
where is the Boltzmann constant. From Eq.(26), we straightforwardly get 

Now let us focus on the modification imposed upon the formalism of 
classical mechanics by the above probabilistic description of stochastic me- 
chanics. 



7 Euler-Lagrange equations 



Let us look at, not the most probable paths (geodesies) of the dynamics 
from A to B which satisfy Eqs.(3) and (4), but all the other paths whose 
unperturbed action is not a stationary. When a system travels along these 
paths, 5A ab (k) 7^ 0. It is either positive or negative. On the other hand, we 
have [6] 



SA ab (k) 



d dL k (t) dL k (t) 



dt dx 



dx 



edt 



(28) 
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where e is an arbitrary variation of x (e is zero at a and b). If 5A a f,(k) > 
(or < 0), we get 



>»<-<«>. 

which can be proved as follows. Suppose J a f(t)edt > and f(t) < c < 
during a small period of time At somewhere between a and b. Since e is 
arbitrary, let it be zero outside At and a positive constant within At. We 
clearly have J a f(t)edt < ce < which contradicts our starting assumption. 
This proves Eq.(29). 

Legendre transformation H = Px — implies [6] 

Fit 

P = (30) 
ox 

Put this relationship back into Eq.(29), we get 

P > (or <) ^ (31) 

for dA a b(k) > (or < 0). In view of the fact that the deviation of P from the 
regular force ^ L is due to the dynamic irregularity, by introducing a "force" 
R of the random perturbation, Eq.(31) can be recast into 

P=^ + R. (32) 
According to Eq.(31), R is positive if dA a b(k) > and negative if dA a b(k) < 0. 



8 Hamiltonian equations 



The total differential of H is 



an f)tf f)fJ 

dH = ^-dP + °-^-dx + -^-dt. (33) 
oP ox Ot 



From Legendre transformation, we can also write [6] 



dH = xdP - ^dx - ^dt. (34) 
dx dt v 1 



11 



Comparing Eq.(33) to Eq.(34), we get 



• 9H , • . ^ dH 

X = ~dP ( ° r <) ~ ~dx } ( > 

which are the Hamiltonian equations for irregular dynamic process following 
the paths whose action is not at stationary. The second equation above 
violates the Newton's second law. Remember that > (or <) is for the case 
where dA ab (k) = ^g^dx - ^^dx > (or <) 0. 

Considering the random force R, Eq.(35) can be written as 

dH 

p> = -is- + * (36) 

This is a Langevin equation. 



9 Energy conservation 

From Eq.(33), we see that if the Hamiltonian of the system is not an explicit 
function of time (^ = 0), we get 

dH dH. dH. 

-dT=dP P + ^ x - (37) 

For the geodesies, using Eq.(3), we get ^ = 0. So the energy is conserved 
during the time evolution of the system. But if the system travels along 
other paths than the unperturbed geodesies, Eq.(35) occurs, we should write 
^jr > (or <) for the paths having dA a b(k) > (or < 0). The energy is no 
more conserved. 

The fluctuation of energy is defined by a 2 H = ((if* — (Hk)) 2 ) = (H%) — 
(Hk) 2 which we will calculate for a segment i of the paths in the time interval 
t = ti — ti-i (see section 5). The action on this segment of the path k can 
be written as 

M = {P iXi - Hi) k r (38) 
So the transition probability is given by 

Pk,i = ^exp[-r]T(PiXi - Hi) k ] (39) 
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where Z{ = J2k ex P [~ r l T (Piii ~ Hi) k ]. After some mathematics, we get 

°h = ^Ys^i-WfzM-wiPiXi-H^] (40) 
z * k 

= -^- + (PxHi)-{Px){Hi). 
r or) 

Let's see an example with free particle. We have in this case (PxHi) = 
2{Hf) and (Pi) {Hi) = 2(H) 2 , leading to 

4 = (41) 

H r dr] v J 



For this kind of systems, (Hi) can be calculated by using Eq.(15) with Fi = 0. 

r_ an {Hi) = mfgl = 1 



We have (Ax, 2 ) = an = = ^ = k B T. Finally, we get 



-I = ^ = 2 ^ ( 42 ) 

Don't forget that we are addressing one dimensional system which can have 
only one particle under thermal noise. For iV-particle systems, the result is 
different. From Eqs.(27) and (41), one can write 



a\ = -2Qk B T 2 (43) 

where Ci = d ^" ! is the heat capacity of the system on the segment i of the 
paths. 



10 Liouville theorem 

Now we look at the time change of state density p(T) in phase space. The 
time evolution neither creates nor destroys state points, so the number of the 
phase points is conserved. The conservation law in a 2-dimensional phase 
space is 

dp d{xpi d(Pp)_ = 

dt dx OP K J 
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which means 



dp dp dp dp ■ dx dP 

Tt = Tt + Yx x+ dP P = - { di + dP )p - (45) 

For the least action paths which satisfy Eqs.(3), the right hand side of the 
above equation is zero, leading to the Liouville theorem 

dJ ^r = °- (46) 

i.e., the state density in phase space is a constant of motion. 

For a system travelling along a perturbed paths having dA a i(k) ^ 0, 
Eqs.(3) vanishes and so does Eq.(46). Considering Eqs.(36) and (45), we 
find 

dp{T,t) dR 

~dT = ~dP p - (47) 

Nevertheless, the determination of the sign of |p is not evident in general. 
In any case, this result means that, in general, if a system does not travel 
along the paths of stationary unperturbed action due to random or molecular 
noise, then its state density may change with time. 

The first consequence of Eq.(47) is that, in general, the phase volume ft 
is no more conserved for irregular dynamic systems since 



dVL d r ,„ , m r dpCT, t) m dR . 

where the average |p is over all the phase volume occupied by the system at 
time t. 



11 Averaged formalism of stochastic mechan- 
ics 

11.1 Euler-Lagrange equations 

It is mentioned above that the dynamics described by the distribution Eq.(10) 
has a stationary average action in association with the stationary path infor- 
mation. In what follows, we use this fact to derive the averaged version of 
the equations of motion of classical mechanics. 
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5A ab can be calculated as follows: 



5A ab = J2Pk( b \ a ) 5A ab(k)+J2 5 Pk(b\a)^ab(k). 
The first term on the right hand side of this equation is 

~ddL k (t) dL k (t) 



Y,Pk{b\a)5A ab {k) = Y.Pk{b\a) f 

J a 

' d dL k (t) 



dt dx dx 
ldL k (t) 



edt 



dt dx 



dx 



edt 



where (■) is the average over all the possible paths. 

The second term on the right hand side of Eq.(49) is 



^5p k (b\a)A ab (k) 



J2^Pk(b\a)\ia[Zp k (b\a)} 



1 

V 

- 5 J2 p k (b\a) In p k (b\a) 



61 



ah 



n 



where we have used J26p k (b\a) = 0. Considering Eq.(12), we get 



' d dL k {t) 
dt dx 



dL k {t) 
dx 



edt = SA ab 



51 



ab 



0. 



V 



Now considering that e is arbitrarily chosen, we obtain 



' d dL k {t) 
dt dx 



dL k (t) 
dx 



0. 



This is the averaged Euler- Lagrange equation. 
From Eq.(30), we can write 



P 



' dL k (t) 
dx 



where 



' dL k (t) 
. dx 



is an averaged force. 
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11.2 Averaged Hamiltonian equations 

By means of Eq.(35) and ^§j— = > the following averaged equations 

exist 



<M§) ^>=-(f )- (56) 

This implies that the mean of the "random force" R over all possible paths 
must vanish, i.e., (R) = 0. 

11.3 Energy conservation law 

From the energy conservation problem discussed above, we see that the only 
Hamiltonian which is physically significant for randomly perturbed hamilto- 
nian systems is the average Hamiltonian (H). Its derivative with respect to 
time (change rate of average energy) is 

The first term of this equation is just 

,dH. ,dH ■ ,dH .. ,dW R , 

<*-> = { 3P P) + = {xR) = <-*-> (58) 

where Wr is a random work done by the "random force" R. 

With the help of the least action distribution pk,% and the definition of 
path information, the second term of Eq. (57) can be given by 

^dj^ H d_p^ p ._^dk 

Y dt Y dt TTj dt K ' 

where Ii = —J^kPk^^Pk^ is the path information on the segment % and in 
the time interval r. Simple calculation leads to 

v (A t A) (60) 



dh 


'd{Ai) 


IdAX 


Tt =r] 


dt 


\ dt 1. 
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On the other hand, 



dp k4 . d(Px) x ^ d{Px) ( . 

> T^PX = : > Pki : 61 

Y dt dt Y dt 

d(Px) _ r dH n . 

d<Pi) . 



Hence we get 



£— = -^ + 2(— x)-2(^)--(A t A). (62) 



This means 

d(H) d(Px) ,dH . . „. 1, A „ N 

- 2^-2<^>-<**>-i(/M) 

etc at r 

where T is the kinetic energy of the system, Wh is the work done by the 
conservative force. We suppose that our hamiltonian system neither gains 
nor loses energy through the irregular random process. So the average energy 
must be a constant of motion, i.e., ^§-- = 0. This means the random work 



must satisfy the following relation: 



This condition seems reasonable because a random force, though statistically 
vanishing, can made non vanishing average work which may change the aver- 
age energy of a system. So if we require that the system statistically remain 
hamiltonian, this work or its power (xR) should satisfy some condition. If 
ever we have (xR) = 0, we can write (^) = 0, meaning that the time 
variation of energy is statistically null. But the average energy may always 
change in time due to the deviation from the least action paths. We have 
the following condition for statistical energy conservation: 



17 



12 Passage from stochastic mechanics to reg- 
ular one 



Of course, when the randomness diminishes, this stochastic mechanics for- 
malism should recover the regular classical one. In fact, the dispersion of 
action or the width of the least action distribution can be measured by the 
variance a 2 = A 2 — A 2 = A 2 — A 2 ab . From Eqs.(5), (8), (10) and (11), we get 

a = —&r> (66) 



A ab = ~hiZ, (67) 

and 

I ab = \nZ + V A ab = \nZ- V ^-\nZ (68) 

When the system becomes less and less irregular, a 2 diminishes and the paths 
become closer and closer to the least action ones having stationary action 
A s ^ b at . When a 2 — » 0, the significant contribution to the partition function Z 
comes from the geodesies having A^ b at , i.e., Z —> exp[—r]A^ b at ]. From Eq.(67), 
A a b — > A s * h at . Then considering Eq.(68), it is clear that the path information 



stat 
ab 



lab — * 0. In this case, the stationary average action Eq.(12) becomes 8 A 
0, the usual action principle, and Eqs.(29) and (35) will recover Eqs.(3) 
and (4). At the same time, diffusion phenomena completely vanish and the 
diffusion laws are replaced by the laws of regular mechanics. The energy 
conservation law is also recovered because Eq.(65) becomes = due 

toA t A = ^-(^) = 0. 



13 Concluding remarks 

The objective of this paper is to present certain mathematical aspects of 
an information method formulated in the attempt to describe diffusion phe- 
nomena within a formalism of stochastic mechanics. We address hamiltonian 
systems undergoing irregular dynamic process due to random perturbation 
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of thermal (molecular) noise and chaotic instability. The dynamic random- 
ness or irregularity is taken into account via the uncertainty associated with 
different paths between two points in phase space and with the different fi- 
nal points from a given initial point. These uncertainties are measured by 
path informations maximized in connection with average action. It results 
that the probability of the irregular process depends exponentially on ac- 
tion (least action distribution). The usefulness of this information method is 
demonstrated by the derivation of diffusion laws such as the Fokker-Planck 
equation, Fick's laws, Ohm's law[10] and the Fourier's law[ll]. 

This information approach with the least action distribution of paths 
implies a deviation from the classical regular dynamics and a formalism of 
stochastic mechanics. The instantaneous behaviors of irregular hamiltonian 
system may violate classical mechanics and its fundamental equations if only 
the unperturbed Hamiltonian is introduced. This is understandable because, 
due to the random noise and instability, exact mathematical treatment of the 
perturbed Hamiltonian is inconceivable even when the system remain statis- 
tically hamiltonian. In this case, only the average Hamiltonian is physically 
significant. 

Another point is that this approach underlines a tight connection between 
the diffusion laws and the maximum dynamic uncertainty associated with the 
average action. Remember that the action here is that of the unperturbed 
regular system, not the real action (does it exist?) of the randomly perturbed 
system. This association of dynamic uncertainty with average action is in fact 
an application of the action principle to the whole dynamic situation, i.e., all 
the paths, not only the most probable paths or geodesies. As a consequence, 
it is the average action, instead of action, which has a stationary under the 
constraint of uncertainty. 

This averaged action principle underlies a formalism of stochastic mechan- 
ics, which can be derived directly from the violated fundamental equations 
by probabilistic consideration. The mathematical form of classical mechan- 
ics can be recovered. The fundamental equations are statistically satisfied. 
It is shown that if the work of "random force" satisfy some condition like 
Eq.(49), the system remains Hamiltonian with a constant average energy. 
We would like to emphasize here that the Liouville theorem is also violated 
by the ensemble of paths or by the instantaneous behaviors of random hamil- 
tonian system. But it seems not obvious to recover this theorem within the 
stochastic version of classical mechanics. This result may have impacts on 
some properties of mechanical systems, e.g., on the Poincare recurrence theo- 
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rem and the mechanical interpretation of the second law of thermodynamics. 
Further work is in progress to tackle this topic. 
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